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We model a melt of monodisperse side-chain liquid crystalline polymers as a melt of comb copoly- 
mers in which the side groups are rod-coil diblock copolymers. We consider both excluded volume 
and Maier-Saupe interactions. The first acts among any pair of segments while the latter acts only 
between rods. Using a free energy functional calculated from this microscopic model, we study the 
spinodal stability of the isotropic phase against density and orientational fiuctuations. The phase 
diagram obtained in this way predicts nematic and smectic instabilities as well as the existence of 
microphases or phases with modulated wave vector but without nematic ordering. Such microphases 
are the result of the competition between the incompatibility among the blocks and the connectivity 
constrains imposed by the spacer and the backbone. Also the effects of the polymerization degree 
and structural conformation of the monomeric units on the phase behavior of the side-chain liquid 
crystalline polymers are studied. 



I. INTRODUCTION 

Side-chain liquid crystalline polymers (SCLCP) are 
comb copolymers composed of a long main chain back- 
bone with elongated rigid side chains regularly attached 
through flexible spacers chain (schematic representation 
on figure nj. Despite their high viscosity, which can be 
two orders of magnitude higher than monomeric liquid 
crystals, SCLCP are good candidates to storage infor- 
mation by locking oriented structures into glassy states 

Q. 

These interesting complexes can be though of as 
supramolecular polymers of rod-coil diblock copolymers. 
In rod-coil diblock copolymers the incompatibility be- 
tween blocks competes against their connectivity lead- 
ing to a frustration of the macroscopic phase separa- 
tion. This competition, together with the tendency of 
the rod segments to align with each other, results in a 
remarkably rich variety of microphase geometries 0, 1^- 
Such patterns are characterized by alternating regions 
of each block with different collective orientations of the 
rods[i|. The main difference between the rod-coil diblock 
copolymers and the SCLCP is the introduction of ad- 
ditional correlations among the rigid molecules brought 
by their linking through the backboneQQ. Based on 
the accumulated theoretical and experimental knowledge 
on rod-coil diblock copolymers 001 nil and comb-like 
copolymers [13 [iH [13 it is natural to expect that mi- 
crophase segregation will also play an important role in 
SCLCP. 

The current experimental and theoretical understand- 
ing of the physical properties of liquid crystalline poly- 
mers derives from the knowledge on monomeric liquid 
crystals and their solutions with polymers. In a pioneer- 
ing work, Vasilenko, Shibaev and Khokhlov devel- 
oped a model for the SCLCP based on Matheson and 
Flory lattice theory for chains with rodlike sections. 
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Figure 1: Model for the side chain polymer. The backbone 
and the spacer are considered as Gaussian chains of rib and Us 
segments respectively while the rigid rods have Ur segments. 
There are nt side groups separated by a distance 5t along the 
backbone. 



The emergence of thermodynamic phases with complex 
nematic orderings was then described as a purely entropic 
effect driven by the main chain fiexibility. This model 
has been later improved by Auriemma, Corradini and 
Va,catello(Tfil| to better include steric repulsion effects. 
These theories neglect the Maier-Saupe anisotropic in- 
teraction and should be thus valid in cases where Maier- 
Saupe interaction between segments is negligible when 
compared to steric repulsion. The study of anisotropic 
interactions in SCLCP was first introduced by Wang 
and Warner Their theory describes the interplay 

between side-chains and backbone nematic orderings by 
means of an extension of the Maier-Saupe theory. In this 
case, the SCLCP melt is treated as a mixture of worm- like 
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molecules and rigid rods with their own nematic order- 
ing tendencies analogously to the model studied by ten 
Bosch, Maissa and Sixou for a solution of nematic 
polymers in liquid crystalline solvents. The Maier-Saupe 
coupling (MS) between the backbone and the side-chains 
is introduced through a constant that describes the effec- 
tive flexibility of the spacer and determines up to which 
extent the side and main chains tend to be perpendicu- 
lar. While in the case of a solution of nematic polymers in 
liquid crystal solvents the MS coupHng between the back- 
bone and the side chains is always positive, i.e., tends to 
align them, in the SCLCP the coupling due to the ef- 
fective spacer attachment competes against the nematic 
coupHng between the side-chain and the backbone. This 
yields an effective nematic coupling between the back- 
bone and the side chains which can be either positive 
or negative, depending on the relative strength of the 
two terms. In Wang and Warner theory excluded vol- 
ume interactions are not taken into account. Thus, we 
should expect it to be valid in the limit where the Maier- 
Saupe interaction is much stronger than excluded volume 
interactions. The possibility of microphase separation 
driven from the competition between the mesogens and 
the backbone + spacer complex is not considered in the 
aforementioned theories. 

In this paper, along the Hues of Holyst and S chick [T^ 
and Reenders and ten BrinkeQ, we study the interplay 
between microphase separation and nematic ordering in 
SCLCP. We treat an intermediate region where both ex- 
cluded volume effects and Maier-Saupe interaction are 
present. We consider that the Maier-Saupe interaction 
acts only between rods and that the excluded volume ef- 
fect acts among all monomers. This should be a reason- 
able approximation when both spacer and backbone are 
flexible. We model the SCLCP as a comb polymer com- 
posed of a backbone and spacer as in reference with 
rigid rods attached to the tips of the comb teeth. The 
comb polymer melt is treated as a set of Gaussian chain 
with Flory-Huggins interaction and Maier-Saupe inter- 
action between rods. The Landau free energy functional 
is derived and the spinodal stability curves are obtained 
as function of the SCLCP structural and interaction pa- 
rameters. Within our formulation the main difference 
between a rod-coil polymer formed by the spacer and 
rod and the SCLCP arises from the inter-rod and inter- 
spacer correlations induced by the connectivity through 
the backbone. Since our main aim in this paper is to 
study the microphase separation instabilities we will not 
consider here the nematic ordering of the backbone. 



II. THE MODEL 

We model a melt of monodisperse Up side-chain poly- 
mers with TV segments of statistical length 6 as a chain of 
rit rod-coil copolymers regularly attached to a backbone 
as schematically represented on figure The backbone 
and the spacer (i.e., the coil part of the rod-coil copoly- 



mer) are flexible Gaussian chains of ni, and monomers 
respectively. The rigid rods are made of rir segments 
linearly arranged along a fixed direction. It is useful to 
define the volume fraction of the rods, fr = ^^^7^^, the 
volume fraction of the spacers fs = ^^^^ff^ and the volume 
fraction of the backbone segments , ft ~ j^. 

There are two important interactions that we consider 
in our model. The first one is the steric repulsion between 
any monomers. The second interaction is the Maier- 
Saupe potential^^ji^l that tends to align the mesogens. 
We write the interaction Hamiltonian as 

(3 Hi = ^ / dxe^/f/^K (x) (x) 



-^Jd^Qf^^ (x) 



Q"" (x) , 



(1) 



where K and K' represent the monomer species indexes, 
i.e, backbone segment, spacer segment or rod segment, 
while exK' and uj are the strength of the interaction 
potentials. Equation Q is already written in terms 
of the volume fractions ipK (x) and the nematic tensor 
Qf^" (x) of the rods, where the local volume fractions of 
the species are defined as 
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where vq is the volume of each segment and j labels the 
side-group position along the backbone. For the back- 
bone, X^ (r) is the curve that describes the conforma- 
tion of the chain a as a function of r, which labels the 
monomers in the chain. For the spacers, separated by 
St = ^ segments, X^'-' (r) describes their conforma- 
tion starting at position X^ ((j ^ ^) <5f:) on the back- 
bone. The unit vector defines the orientation of the 
rod j, connected to the end of the spacer at position 
X^.^(n,) + X^(0--i)Jt). 

The nematic tensor of the rods is defined as 
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For the sake of simplicity we will considered that the 
backbone and the spacer are made of the same segments. 
In this way, they can described by a single field (fic = 

<i>b + <^s- 

Following Gupta and Edwards [IJl we rewrite the par- 
tition function of the entire melt in terms of collective 
variables 4>c,4>r,Q as 

Z = jT^<f>cj Vcf>r j PQe-'^^^[*-'^'-''^lZo[(/)e>r,Q] 

(6) 

where Zq is the partition function of independent SCLCP 
molecules which yields the entropic contribution to the 
total free energy and is given by 

x,5(0,-0,)<5(q-q)5(K| -1) 

^^-PHo[K,^-^] . (7) 

The Hamiltonian of unperturbed Gaussian chains, Hq, 
is written as 



2Nb^ 
Vt 



Vt 



2' 



dr 

As discussed by Shinozaki et al.^3|) this assumption 
of Gaussian statistics for the chains is probably not the 
more reaHstic since the chains might be stretched rela- 
tive to the Gaussian case. Nevertheless, for comb with 
an evenly spaced teeth and small monomer density, the 
Gaussian description is still reasonable. 

We now use auxiliary fields he, hr and hg" to express 
the Dirac deltas in equation Q, 

-i f drhr<t>r -i f drh^"'Q^"' 



exp{-Fi [hr,hchQ]} 



(8) 



where 
exp{-Fi [hr, hchg]} 



TT"" tt"* 



Following references |2l| and (23|, the integrals over 
the auxiliary fields in equation © are approximated ac- 
cording to the steepest descent method. We write the 
entropic contribution to the free energy, Fi [hr,hchQ], 
as : 



Fi [hr, hchQ 



q{2) pu per 
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where G^'^^^'^'"^ is the matrix of non interacting pair cor- 
relation functions : 



(11) 



The non-interacting chains pair correlation functions 
are given by 



Gkk- (x,x') = (0k(x) ^k' (x) 
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where (...)q means the average over the partition function 
of the unperturbed Gaussian chains |S1 and (I>k (x) and 
Q'^'^ (x) are defined in equations JSJ , (EJ , © and 101 . 
We can now expand the total free energy of the system 
as a function of the order parameters^ and Q following 
the same procedure as in references j22| and Q- At this 
point, we assume that the system is incompressible, i.e.. 



(12) 



(j)c (x) -t- 0r (x) = 1 , 

which for q 7^ means 

After applying condition in Hamiltonian, the excluded 
volume interaction term reduces to 



X J dx(j}{x) , 

where x =j_%c + £rr — '^^cr is the usual Flory- Huggins 
parameterj23. 

The second order expansion term, which defines the 
stability limits of the isotropic phase against fluctuations 
of the order parameters, is given by 



(13) 
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The second order vertex r^^^'"^'"^ is the inverse of the 
matrix G^^^"^ , as calculated in references and plus 
a diagonal matrix that contains the terms for the Flory- 
Huggins interaction and the Maier-Saupe potential: 
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This approximation for the second order vertex, known 
as random phase approximation (RPA) is considered as 
reasonable for dense mixtures of strongly correlated poly- 
mer chains, where the fluctuations are very small (2^|2fl|. 
The calculation the matrix Ijl4|l follows along the same 
lines of reference Q- The major difference here is that 
there are additional correlations between the rods and 
the spacers attached to the same backbone, in contrast 
to the rod-coil copolymer studied in reference j3| • For in- 
stance, the Fourier transform of the rod-rod correlation 
function in the SCLCP is given by 
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where, Krr {x) = ^ [cos {x) — 1 + x Si (x)], is the rod- 
rod correlation function for the rod-coil copolymer stud- 
ied in reference 0, Fr {x) = ^^^^ is the form factor of a 

_ 2 / 2 \ 

e +sinh(^)n — 1 



rod of length I and D„ {x) 



is the 



Debye function for the backbone. Here we use I = rirb as 

the length of the rod and Rs = -\/^ and Rb = -y/^ 
as the gyration radius of the spacer and of the backbone 
respectively. 

Note that in the limit rit = I and St the correlation 
function Ijl5|l reduces to the rod-rod correlation function 
for the rod-coil copolymer studied in reference 01 . In the 
opposite Hmit, for nt ^ 1, which is the case of SCLCP the 
rod-rod correlation function is dominated by the inter- 
rod correlations within a SCLCP and it is mostly influ- 
enced by the backbone and spacer correlations. This is 
what leads to the main aspect of the physical proper- 
ties of SCLCP when compared to rod-coil copolymers 
and to the strong dependence of the mesophases on the 
structural parameters of the SCLCP (a|0|- The other 
correlations functions are shown in appendix El 

Also, as done in reference 0, we consider a nematic 
ordering parallel to the wave vector q: 



Q""" (q) = Q (q) 



qfq" 



(16) 
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Figure 2: Effect of the polymerization degree and the rela- 
tion between the Maier-Saupe and Flory-Huggins interactions 
on the phase diagram for SCLCP for (a) ^ = 1 (b) — = 8. 
The isotropic phase is denoted by "I", while the microphases, 
nematic and smectic phases are labeled as "M", "N" and "S" 
respectively. The diagrams were calculated for ^ = 25, dif- 
ferent values of and nt and St as indicated . 



III. PHASE BEHAVIOR 

A. Stability limits of the isotropic phase: phase 
diagram 

The stability of the isotropic phase against fluctua- 
tions of the order parameters is given by the spinodal 
curve. The roots of the determinant of r(2)A'J'P<^ define 
an expression for x as a function of the wave vector 
the structural parameters n^, rig, nt, St and the ratio 
between Maier-Saupe and Flory Huggins interactions: 



det 



^ ^ x[ q,n.r,ns,nt,St,-j . (17) 
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For rir, rig, rit, 5t and ^ fixed, the spinodal curve is de- 
termined by the minimization of x with respect to q. This 
has been calculated using a program written for Math- 
ematica and available from the authors by request. As 
discussed in reference j23| , if the value q* that minimizes 
equation ifTTIl is equal to zero, the isotropic phase is un- 
stable to a nematic phase. If q* 7^ 0, then the isotropic 
phase is unstable against modulations of the volume frac- 
tion, which gives raise to microphase separation. On that 
basis, it is possible to construct a phase diagram that 
sketches the the transition from an isotropic phase to an 
ordered stated as a function of the structural parameters 
of the model. For shortness and simplicity, in this paper 
we denote as IN the spinodal obtained from the condi- 
tion g* = and IM the spinodal calculated according to 



B. Effect of the polymerization degree 

In figure |2l we show the effect of the polymerization 
degree on the phase diagram for two different values of 
the ratio ^. In both cases, the phase diagram for St = 1, 
nt — 1 is the same as the rod-coil diblock copolymer 
phase diagram calculated in reference 01 ■ The region 
under the curves corresponds to an isotropic phase, de- 
noted as 'T' in the figure. Above them, the phase can 
be either nematic (N), microphase (M) or smectic (S), 
which is superposition of the nematic ordering and den- 
sity modulation. When the Maier-Saupe interaction is 
of the same order of the Flory-Huggins interaction, i.e. 
^ = 1, the phase diagram shows the possibility of having 
microphase segregation before the transition to a smec- 
tic state, like in phase diagram (a) of figure |21 If the 
value of the ratio ^ is such that there is no crossing 
point between the two spinodals, as in phase diagram 
(b) of figure |21 there is no microsegregation for any vol- 
ume fraction. This might be attributed to the fact that 
in this case, the Maier-Saupe interaction is stronger than 
the Flory-Huggins interaction. Therefore therefore the 
nematic ordering occurs first than the modulation of the 
density. 

We note that the spinodal curves for the isotropic- 
nematic instabilities, are given by 



IM — 7 . 1101 

) 

This instability line does not depend on the polymer- 
ization degree of the SCLCP but only on the ratio i.e, 
the total number of monomers in each rod-coil comb plus 
St. This is exactly the same dependence as the case of a 
single rod-coil studied in reference . It is a consequence 
of the approximated model we adopted where the rods 
only interact through the Maier-Saupe coupling among 
themselves. Since the isotropic-nematic instability occurs 
at the wave- vector q = 0, the inter- rod coupling along a 




0.2 



0.3 0.4 0.5 0.6 0.7 0.8 

fr 

Figure 3: Value of the wave vector q* that defines spinodal 
for the formation of microphases.The structural parameters 
were set as nt = 30, St = 3, and rim = 25. 

SCLCP chain is not influenced by the backbone or the 
spacer. 

The spinodal curves for the isotropic-microphase in- 
stabilities are shifted towards higher temperatures as the 
degree of polymerization increases. According to Finkel- 
mann ad Rehage [23|, this behavior is well established 
experimentally, and can be understood by the restriction 
of the translational and rotational motions of the meso- 
genic molecules due to the linkage to the backbone. After 
rit ~ 10, the chain length is such that there is no correla- 
tion between the first monomer in the chain and the last 
added. Therefore, the instability curve is hardly affected 
by the polymerization degree. 

The spinodal curve IM is also infiuenced by the ratio ^ 
, specially in the rod-rich regions, as can be seen compar- 
ing phase diagrams (a) and (b) of figure[2 However, this 
dependence vanishes with increasing nt and after nt = 5 
the spinodal curve IM is not longer affected by the ratio 

X ■ 



C. Characteristic length 

In figure El we show the value of q* that defines the 
spinodal curve IM as a function of the rod volume frac- 
tion. The solid line is the wave vector for the length of 
the rod-coil part of the SCLCP, i.e. qrs = ^^y^ - it 

can bee seen from the figure, even in the rod-rich region, 
the modulation period of the microphase is determined 
by both rod length and spacer conformation. Neverthe- 
less, up to fr w 0.8, q* is always smaller than qrs ■ As 
the polymerization degree is increased (see dotted and 
dashed lines in figure EJ, this difference increases. On 
the basis of this information and the experimental data 
on SCLCP discussed by Noirez and coworkers in refer- 
ences '28j[29j^33|, we suggest the microphase region is 
formed by structures in which the backbone is confined 
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between layers containing the rigid cores. Therefore, the 
difference between q* and qrs might be attributed to the 
conformation of the backbone. 



D. Temperature equivalence 

So far, all the information we have discussed is based 
on the calculation of the Flory-Huggins parameters as a 
function of the structural parameters of the model and 
the ratio ^ . Since it is well known that the Flory-Huggins 
parameter has a dependence on the inverse of tempera- 
ture of the form 



T{K) 



A „ 



(19) 



it is interesting to relate the constants A and B to the 
SCLCP chemical constitution in order to give a more 
realistic characteristic to our results. 

From thermodynamic arguments applied on the 
Flory's lattice model for a polymer in a solvent, it can be 
shown that the temperature dependent term of expres- 
sion is the enthalpic contribution from rearrange- 
ment of the contacts polymer-polymer, solvent-polymer 
and solvent-solvent . See for example reference [23l for a 
review on this subject. This term can be estimated us- 
ing the Hildebrand solubility parameters of the polymer 
and the solvent segments. Following the enthalpic 
contribution to the Flory-Huggins parameter is 



A = 



R 



(^^pol ^solv) 1 



(20) 



where the 5' s are the Hildebrand solubility parameters of 
each component of the mixture, Vseg is the common seg- 
ment volume and R is the ideal gas constant. The Hilde- 
brand solubility parameters are approximated through 
the relationship with the cohesive energy and the group 
contribution method j^J] . 

The athermal term in expression Ijl9|l is related to in- 
ternal degrees of entropy of the polymers. According 
to different authors it might have its origin in features 
such as conformational asymmetries of the components 
of the mixture and nematic interaction between 

the segments originated by the flexibility differences be- 
tween the polymers. In general, this term is determined 

empirically • 

This equivalence between temperature and the x Pa- 
rameter allows us to compare our results with experi- 
mental data available in literature. As an example, in 
appendix El we calculated the Flory-Huggins parameter 
as a function of temperature for the polyacrylates shown 
in figure El with Ug = 2, rit = 30, St = 2 and = 9. 

This SCLCP has a isotropic-nematic transition at 
Tj^^ = A25K m. As a test for the validity of the this 
correspondence between temperature and x parameter. 
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Figure 4: .Transition temperatures as a function of the spacer 
length. Parameters:(5f = 3, Ur = 10, nt = 30 



we calculated the transition temperature for this poly- 
mer using our model. Setting ^ = 4.6 and B = —0.05 we 

obtain a critical temperature of = A30K. Since we 
have done several approximations in the model, and the 
results depend on two adjustable parameters, we must 
not be too optimistic about the agreement of the model 
presented here with respect to the experimental data. 
Nevertheless, the agreement achieved can be considered 
as a good indicative of the possibility that main physical 
features of the systems are present in our phenomenolog- 
ical model. 

From now on, we will use relation lfT9|l to talk about 
temperature instead of the Flory-Huggins parameter. In 
this way our discussions will have a more intuitive back- 
ground. 



E. Influence of the structural parameters on the 
phase diagram morphology 

According to Finkelmann and Rehage j the 
schematic phase behavior of SCLCP follows the same 
trends as low molar mass liquid crystals. For a fixed 
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Figure 7: Transition temperatures as a function of the meso- 
gen length. The parameters were set as St = 3, ris — 10, 
nt = 30. 



phase is not necessarily favored by increasing . In fig- 
ure El we show the phase diagrams for a larger value of 
St. As it can be deduced from the figure, the behaviors 
described above are maintained despite the larger value 
of St. 

As discussed by Auriemma|Tfil| et al. and by Finkel- 
mann and Rehage j23|, another characteristic behavior 
of low molecular weight liquid crystals and enhanced in 
SCLCP, is the increase of the isotropization temperature 
with increasing length of the mesogenic groups. This 
trend is also present in our model. In figure , we show 
the transition temperature as a function of rir for differ- 
ent values of the ratio ^. In this case, the smectic phase 
is favoured by both the increase of rir the interaction 
ratio ^. 



Figure 6: Polyacrylate used as example to calculate the de- 
pendence of X on temperature. 



mesogen length, increasing the length of the fiexible 
spacer or the length of the backbone segment neighboring 
side chains decreases the instability temperature towards 
nematic or smectic mesophases. Some aspects of our re- 
sults are however different. First, in our model, for a 
fixed St, the smectic phase is favored as is increased 
only if w < tOc- Here, lUc is a function of the isotropic- 
microphase transition temperature curve. It is defined 
as the curve Wc that crosses the curve IM at the point at 
which it starts decreasing with increasing Us. For larger 
Maier-Saupe couplings the general trend is opposite to 
that related by Finkelmann and Rehage. In figure we 
show the phase diagram as a function of the spacer length 
for two different values of — . In the first case, the behav- 
ior is as described by Finkelmann and Rehage, while in 
the second case, where the ^ ratio is larger, the smectic 



F. Microphase segregation 

As mentioned before, an important feature of our 
model, is the possibility of microphase segregation be- 
sides the formation of nematic and smectic phases. In 
this sense, our phase diagrams for SCLCP are more sim- 
ilar to rod-coil diblock copolymers phase diagrams. In 
these systems, the mutual repulsion of the blocks, due 
to the difference in chain rigidity, and the constrains im- 
posed by the connectivity of the blocks result into the 
formation of supramolecular structures as small as few 
nanometers0- If we think about the SCLCP as being 
a polymer in which the monomeric units are rod-coil di- 
block copolymers, it is quite natural to expect the exis- 
tence of such microphases. 

This kind of interplay between liquid crystalline 
mesophases and microstructure domains is documented 
experimentally for SCLCP made of diblock copolymers, 
as in references js^H^I . In these SCLCP, the repeated 
unit is a homogeneous diblock copolymer with a liquid- 
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crystalline monomer grafted to one of the blocks. Such 
kind of SCLCP exhibit microdomains with spherical, 
hexagonally packed cylindrical and lamellar structure. 
The morphology of their phase diagram can be related 
to the molecular weight, the block composition and the 
X parameter. 

Therefore, the existence of microphases in SCLCP is 
a possibility that should be considered when performing 
experiments and analyzing experimental data. 

It is worth to note that despite our model predicts mi- 
crophase segregation for SCLCP with certain structural 
characteristics, we can not differentiate the symmetries 
of this microphases. To distinguish between the different 
possible symmetries, it would bee necessary to perform 
calculations with higher order terms of the free energy, 
as done by Reenders et al. Q for the case of rod-coil 
diblock copolymers. 



G. Frank elastic constants 

As explained in the review written by Stephen and 
Straleyj2y|, the gradient terms in Q of the free energy 
are related to the Frank elastic constants. 

The usual expression for the Frank elastic energy in 
terms of the director n is 



Feiasuc = ii^i (V • n)' + ii^2 (n • V X n)' 

+ iA-3 (n X V X n)' . (21) 

The constants Ki are the splay, twist and bend Frank 
elastic constants. In general, the bending constant is 
larger than the other two, which are about the same order 
of magnitude. In many cases equation lf2T]l is too complex 
to be of practical use. Thus, it might be useful consider 
all the constants as equal. In that case the elastic energy 
is reduced to 



Equation ^22\ is not quantitatively correct but gives a 
qualitative idea of the distortions in nematic Hquid crys- 
tals 113. 

In figure IHl we show the correlation function for the ne- 
matic tensor calculated near the isotropic-nematic tran- 
sition temperature. This correlation function can be ap- 
proximated by a function of the kind 



Gqq (9) 



(23) 



The inverse of Gqq corresponds to the expansion of 
the free energy term proportional to the square of the ne- 
matic order parameter in power series of the wave vector. 

Gqq [q) x 10^ 
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Figure 8: Correlation function for the nematic tensor calcu- 
lated for 5t — 2,, Us = 7,nr = 10, nt = 30, ^ = 4 and 

T = l.lTcrit- 



FeiasUc ~ -K diUjdjrii, 



(22) 



The coefficient ^ carries the dependence on the Maier- 
Saupe coupling. When /i = 0, the correlation function 
diverges at g = 0. Therefore, is defined by the the IN 
curve given by equation ltT8|l . 

_ 5 2uj N 



where the sum over repeated indexes is employed. 



The coefficient L , which corresponds to 



nt 


nl (137-f 28nt) 


[ris + 5tf + 220ns (n^ + in^St 


+ 3St^j 


7 


3ri,2 [us -t- 


5tf + {nl + Sn^St + 35t^) 





(24) 



depends only on the structural parameters nr,ns,St, 
and nt . Regarding to the form of Ij23|l , we can interpret 
the quantity = ^ as a correlation length related to the 



nematic ordering. Also, coefficient L is the gradient term 
of the free elastic energy and thus can be associated to 
the elastic constant K. Note that this elasticity is purely 
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K{nt) 
K(l) 
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nt 

Figure 9: Elastic constant K as a function of the poly- 
merization degree for various values of ris at a temperature 
T = l.lTcrit- The other parameters were set as fr = 0.4, 
(5f = 3, ^ = 4, nt = 30. 
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Figure 10: Volume fraction correlation function for tempera- 
ture close to the transition. The parameters used are St = 3, 
Ur = 10, Us = 12, and ^ = 4 . 



entropic. Therefore, by studying the behavior of the cor- 
relation function Gqq with the degree of polymerization 
and the structural parameter of the SCLCP, its possible 
to obtain information of how the polymerization degree 
and the size of the spacer affect the nematic ordering. 

In figure El we show the coefficient K as a function of 
the polymerization degree for various values of rig near 
the isotropic-nematic transition temperature. Clearly, 
the additional correlations introduced polymerization de- 
gree raise the energy cost of the nematic distortions. This 
reffects as a rise of the nematic correlation length. From 
ll2^ it is clear that K has an asymptotic behavior of 



nf. Therefore the nematic correlation length scales as 
nt. While in low molar mass liquid the elastic energy 
represents a small fraction of the overall energy j^S], its 
contribution can be amplified considerably by the poly- 
merization in SCLCP . 

On the other hand, the increase of the spacer size 
slightly reduces the value of K. The longer the spacer, 
the weaker the coupling between the mesogens and back- 
bone conformation and the shortest the nematic correla- 
tion length. 



H. Correlation length for the density fluctuations 

Information about the correlation length can be ex- 
tracted from the volume fraction correlation function in 
the reciprocal space, G^^ (g). This function is calculated 
as the first element of the inverse of matrix r^^^'^"'"^. 
In figure nni we show df,^ {q) for different temperatures. 
Even at very high temperatures, its maximum is always 
at q* ^ 0, i.e., the homogeneous state is not the lowest 
free energy state. Therefore, we suppose that near the 
transition temperature, the correlation function G^^ (q) 
is proportional to a function of the kind 



54,4, (q) 



For temperatures very close to the transition temper- 
ature, we can use the approximation 



G4"l> (9) 



1 



where go is the wave vector that maximizes C/^^ (g) , and 
^ is the correlation length. 

In figure we show the correlation length ^ as a func- 
tion of the reduced temperature t = ^^^^ for different 
degrees of polymerization. The change in the value of the 
correlation length from nt = I to nt = 2 clearly shows the 
the importance of the correlations induced by the linkage 
of the mesogens to the backbone. As the polymerization 
number increases, the correlation length increases. As 
discussed by Finkelmann and Rehage (23 , after a certain 
number of monomers added to the SCLCP, the correla- 
tion between the last and the first monomer is lost. 

As expected for a second order transition, the correla- 
tion length diverges at i = 0. Nevertheless, the form of 
the correlation function is such that all the allowed wave 
vectors form a spherical surface of radius go. Such a con- 
dition is characteristic of systems with a fiuctuation in- 
duced first order transition, as deduced by Brazovskiijs^. 
Higher order calculations would be necessary to corrob- 
orate this statement. 



IV. CONCLUSIONS 

As discussed above, the model we have proposed for 
SCLCP is in good agreement with the experimental ob- 
served trends discussed by Finkelmann and Rehage j23| 
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Figure 11: Correlation length as a function of the reduced 
temperature for different values of the polymerization de- 
gree. The parameters used are St = 3, rir = 10, Us = 12, 
^ = 4 and nt = 30. 

and Shibaev and PlateQ for the transition temperatures 
of the SCLCP. Moreover, for and adequate value of the 
ratio ^, our model predicts the existence of microphases 
with a modulated volume fraction but without orienta- 
tional order . These microphases, originated from the 
mutual repulsion between the rigid and flexible blocks of 
the SCLCP, are already observed in SCLCP made of di- 
block copolymers as explained in references j^ElH^- Nev- 
ertheless, further calculations are necessary in order to 
identify their symmetry . 



As expected, the additional correlations introduced by 
the linkage of the mesogenic groups to a polymeric back- 
bone through the flexible have a strong influence on the 
SCLCP critical behavior. This statement is supported 
by two main features. The first one is the swelling of the 
characteristic length of the microphase domains, that we 
have attributed to the backbone conformation. The sec- 
ond one is the the infiuence of the spacer length Us, and 
polymerization degree rit on the phase diagram topology, 
correlation length of the volume fraction correlations and 
elastic energy associated to the nematic distortions. 
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Appendix A: THE NONINTERACTING 
CORRELATION FUNCTIONS 

The correlation functions are calculated as explained 
in references with the difference that there is a 

rigid rod attached at the end of each tooth of the comb 
polymer. 

As an example, we calculate the correlation function 
rod-rod. This function has two contributions: one that 
comes from the correlation of the segments within the rod 
and the other one coming form the interaction between 
segments of two different rods. In this way: 



— Krr (Nrbq) + e-^^'^^Ffl {Nrbqf (On, (qRb) - —) 

nt \ nt J 



(Al) 



where, Krr{x), Fr (x) and D„ (x) are the the rod- rod 
correlation function, thee form factor of a rod of length 
I and the Debye function for the backbone as explained 
in the text. 

The other correlation function, calculated in the same 
way, are 

Gbb (q) - f^Dn, (qRb) , 

G,,(q) = !lDnAqRs) + 

nt 



f^Fn, {qRsY {qRb) - - 

nt 



where (qR^) 
of the spacer. 



is the structure factor 



Grs (q) = frfsFr (qbur) {qR^ 



— + ( {qRb) 
nt 



1 



nt 
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6 

1 _ eT^^'''^*F„, (qR,) 



Gcc (q) = G,, (q) + (q) + (q) , 



and 



is the structure factor 



"6 1 9 -f^b 

e "b -1 

of the backbone segments between two side groups, 



Gsb (q) = 2cfi,fsF„^ (qRs) 



1 _ e^'^'^'^'Fr,, (qRb) 



where we have defined A^'' = 2-| 



G^". (q) = frfs^^^F^^ {qR,) KQ^r {rirbq) 



X < he 



(qRb) - 



nt 



b'^q^rii, 



1 _ eT^?''''**^,,, 



where Fq__ (a;) = [3 sin (a;) — 3a; cos (x) — x^S'i [x) 

1 j>- / \ 4a: — a:; cos a: — 3 sin X— x^Sifa:) i r» i • r 

and Aqj. (x) = ^ as defined m ref- 
erence Q, 



Gqc (q) — Gqs + GqIj. 
Appendix B: THE HILDEBRAND PARAMETERS 

According to reference for non polar substances, 
the solubiHty parameter can be calculated through its re- 
lationship with the cohesive energy and its molar volume: 



V 



(Bl) 



Using to the group contribution method, the solubility 
parameter of a monomer unit is expressed as the sum of 
independent contributions from the atomic groups that 
constitute the unit. 



Ecoh — ^ 



coh ' 



(B2) 



In the case of interest, we are considering the mixture 
of two components: the comb polymer (backbone and 
spacer). Therefore, it is convenient to choose the the 
segment unit as having a molar volume equivalent to that 
of the CH2. In this way, using ijBljl and the with the data 
in reference is [131, we find that the Hildebrant solubility 
parameter for the fiexible part of the polyacrylate is 



Spot = 16.13 (Jcm^^) 



G^Q^'^(q) = f^Y.'^r 



pa- 



—K 

nt 



(2) 

Si 



+e-2i,R. fp^^ (d^^ (qR,) - 1-) 5,- 



where the tensors T/^'^'P'^ and the functions Ksi are de- 
fined as in reference 

Therefore, the correlations used to calculate matrix ITTl 
are defined as follows: 

Gcr (q) = Grs (q) + Grb (q) , 



For the rigid part, we must divide the cohesion energy 
of the whole group by n^. This is equivalent to substitute 
the rigid rod by a chain of Ur units of GH2 each one with 

a solubility parameter Ssoiv = V ''^ . In this way we 



obtain the Hildebrant solubility parameter for the rigid 

core, 



7.73 (Jem 



-3^l/2 



By substituting and in 6poi and Ssoiv inlJOl it is possible 
to calculated transition temperatures. 
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